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>>■ Abstract 

A basic feature of fluid mechanics concerns the frictionless phase-space dynamics of particles in 
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an incompressible fluid. The issue, besides its theoretical interest in turbulence theory, is important 
in many applications, such as the pollutant dynamics in the atmosphere, a problem relevant for 
civil protection monitoring of air quality. Actually, both the numerical simulation of the ABL 
(atmospheric boundary layer) portion of the atmosphere and that of pollutant dynamics may 
generally require the correct definition of the Lagrangian dynamics which characterizes arbitrary 



fluid elements of incompressible thermofluids. We claim that particularly important for applications 
£-v ', would be to consider these trajectories as phase-space trajectories. This involves, however, the 

unfolding of a fundamental theoretical problem up to now substantially unsolved: namely the 
determination of the exact frictionless dynamics of tracer particles in an incompressible fluid, 
treated either as a deterministic or a turbulent (i.e., stochastic) continuum. In this paper we intend 
to formulate the necessary theoretical framework to construct such a type of description. This is 
based on a phase-space inverse kinetic theory (IKT) approach recently developed for incompressible 
fluids (Ellero et al, 2004-2008). Our claim is that the conditional frictionless dynamics of a tracer 
particles - which corresponds to a prescribed velocity probability density and an arbitrary choice of 
the relevant fluid fields - can be exactly specified. 

PACS numbers: 47.10.ad,05.20.Dd 



I. INTRODUCTION 

This work (together with Refs.pj) is a part of a research project related to the theoretical 
description of pollutant dynamics in the atmosphere, a subject relevant for civil protection 
monitoring of air quality in the atmosphere. Here, we refer in particular to the so-called 
ABL (atmospheric boundary layer) portion of the atmosphere where the earth surface (land 
or water) has a direct influence and most of pollution releases occur. In fact, the ability to 
predict the dynamics of anthropogenic pollutants, especially in order to estimate their con- 
centration at ground level, is a prerequisite for environmental investigations. A critical issue 
is therefore the identification of mathematical models, able to give reliable predictions for 
pollutants concentrations in the presence of complex terrain and for prescribed weather pro- 
files. This involves the ability to simulate pollutant dynamics in a variety of different physical 
conditions. In fact in the ABL, sufficiently close to the ground, the atmosphere is mainly 
characterized by a turbulent flow arising from the wind shear produced by friction with the 
ground surface. Instead, at the top of the ABL, in the free atmosphere, the wind speed 
is approximately geostrophic and therefore possibly laminar. The stronger the wind, the 
more intense is the generated turbulence arising close to the ground, whose properties may 
be very different, depending on the vertical temperature gradient of the atmosphere. Since 
turbulence reinforces mixing it tends to homogenize the fluid much more quickly than would 
a laminar flow, thus preventing local accumulation of pollutants. In addition meteorological 
parameters (of the atmosphere) are strongly affected by the earth's surface through dynam- 
ical processes (friction of the air over the surface and through thermal processes heating or 
cooling of the air in contact with the ground). Until recently, fully reliable mathematical 
methods of this type, able to take into account the full complex phenomenology of the at- 
mosphere and simulate the dynamics of particle or gaseous pollutants in the atmosphere, 
have been missing. For this reason in the past, most of predictions for pollutant transport in 
the atmosphere have been based on wind-tunnel experiments. Purpose of the presentation 
is an overview of the mathematical models currently available and a brief analysis of new 
theoretical developments in the field. In particular, here we intend to refer to the statistical 
approach developed by Marco Tessarotto et al. [lj based on the IKT (see also Ellero et 



al., 2000-2008, [2, 



i, y, y, y, y, 



261] ). In the sequel, we shall concentrate on those issues that 



we consider mostly relevant for the present investigation. We intend to show that the IKT 



approach is a useful theoretical framework that can be used to simulate the dynamics of 
pollutants in prescribed fluid flows, taking into account, besides the velocity and pressure 
fields also the temperature of the same fluid. The mathematical approach appears therefore 
susceptible of actual applications for simulations of pollutants dynamics in the atmosphere, 
such as those relevant to civil protection monitoring. 

II. THE PRESENT STATE OF RESEARCH 

The subject of the investigation (pollutant dynamics) involves some of the most disputed 
and still open issues in fluid dynamics. The motion of particles in a fluid as they are pushed 
along erratic trajectories by fluctuations of the fluid fields is fundamental to transport and 
mixing in turbulence. It is essential, for example, in combustion processes [7j, in the indus- 
trial production of nanoparticles [8j as well as in atmospheric transport, cloud formation 
and air-quality monitoring of the atmosphere [9|, llO]. The Lagrangian approach - denoted 
as Lagrangian turbulence (LT) - has been fruitful in advancing of the understanding of the 
anomalous statistical properties of turbulent flows [ll]. In particular, the dynamics of par- 
ticle trajectories has been used successfully to describe mixing and transport in turbulence 
3, Qjj. Nevertheless, issues of fundamental importance remain unresolved. As it is well 
known, the description of fluids can be performed choosing either Eulerian or Lagrangian 
viewpoints. The two approaches are equivalent and, if fluid dynamics were fully understood, 
one should be able to translate Eulerian properties into Lagrangian ones and viceversa. We 
are still quite far from this goal. 

In the past the treatment of Lagrangian dynamics in turbulence was based on stochastic 



models of various nature, pioneered by the meteorologist L. F. Richardson 13]. These 
models, which are based on tools borrowed from the study of random dynamical systems, 
typically rely - however - on experimental verification rather than on first principles. A 
first example is provided by refined stochastic models - based in non-extensive statistical 
mechanics [l4J, ll5[] - which try to reproduce the observed intermittency of turbulence 16, 
17|. Another approach is based on Brownian- motion-type models, which involve stochastic 
particle equations of motion and an associated Fokker-Planck equation (see, for example, 
Subramanian et al.,2004 }18|). In such models the diffusion tensor in the Fokker-Planck 
equation may depend on the local velocity shear rate, implying a possible violation of the 



classical (equilibrium) fluctuation-dissipation theorem [19j]. A further class of models has 
been based on kinetic theory, enabling us to obtain the Fokker-Planck equation for the 
motion of Brownian particles in rarefied nonuniform gases. In particular, by expanding 
in term of the mass-ratio (of the light-gas molecular mass to the Brownian-particle mass) 
the Boltzmann collision operator, the light-gas distribution function was approximated by 
the first two terms in the Chapman- Enskog expansion. As a result, it was found to obey a 
bkker-Planck equation with a diffusion tensor independent of the light-gas velocity gradients 
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2l|. 



However, in most cases there remains a lack of experimental data to verify the reliability 
of such models. Verification can be based, in particular, on the measurement of fluid particle 
trajectories, obtained by seeding a turbulent flow with a small number of tracer particles and 
following their motions with an imaging system. On the other hand, the accurate evaluation 
of the Lagrangian velocity in laboratory turbulence experiments requires measurements of 
positions of tracer particle by using a suitable tracking system able to resolve very short time 
(and spatial) scales. In practice this can be a very challenging task since particle motions 
must be measured on very short time scales of the order of the Kolmogorov time, 

r = (u/e)^ (1) 

where v is the kinematic viscosity and e the turbulent energy dissipation. Similarly, to get the 
Lagrangian acceleration one should have experimental access to time scales again comparable 
or even smaller than the Kolmogorov time scale of the flow. Another important physical 
observable is the pressure (and the related pressure gradient), which is typically hard to 
measure experimentally. Despite these difficulties, growing interest in studying Lagrangian 
turbulence is motivated by the recent advances in laboratory and numerical experiments 
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24j). In particular, the role of numerical simulations represents a new challenge, 
since they are expected to represent alternative tools to laboratory experiments. 

As for the theory itself, rigorous results have been scanty, probably because of the subject 
complexity. In particular, an open issue is the very definition of the dynamics of tracer 
particles which may injected in the fluid with arbitrary initial velocities. These velocities in 
practice may be generally very different from the local fluid velocity. In contrast, in customary 
treatments the dynamics of tracer particles is identified with that of the fluid particles, which 
is determined uniquely by the Navier-Stokes equations. As a consequence, the dynamics 
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(i.e., in particular the acceleration) of tracer particles in incompressible thermofTuids might 
be very different from that of the fluid elements. This problem is not only a challenging 
theoretical issue, but is obviously of primary importance for its implications in computational 
and environmental fluid dynamics. 

In this connection, however, recently there has been an important breakthrough, repre- 
sented by the discovery of the so-called inverse kinetic theory (IKT) approach for incom- 

|25| ] the same approach has been 



pressible fluids J2J, y|, |J, |5|, |6( . In accompanying papers |l|, 
extended to the treatment of incompressible thermofluids. 

In this paper we claim that the IKT of Ref.]2J permits us to determine uniquely the 
(frictionless) dynamics of tracer particles. 

In fact, a key feature of the this type of approach is that it affords a straightforward con- 
nection between Eulerian and Lagrangian fluid descriptions. This is achieved by identifying 
the relevant fluid fields, which are assumed to be defined in a suitable domain Q C R 3 (fluid 
domain), with appropriate moments of a suitably-defined kinetic distribution function den- 
sity /(x,t) [with x = (r, v) 6 T, x and T denoting a suitable state-vector and an appropriate 
phase-space] which is assumed to advance in time by means of Vlasov-type kinetic equation. 
In such a case, the time-evolution of the kinetic distribution function is determined by a 
kinetic equation which, written in the Eulerian form, reads 

L/(x,t) = 0. (2) 

Here L is the streaming operator Lf = £/ + J^ • {X(x, £)/} and X(x, t) = {v,F(x, t)} 
a suitably smooth vector field, while v and F(x, t) denote respectively appropriate veloc- 
ity and acceleration fields. As a main consequence the approach can in principle be used 
to determine in a rigorous way the Lagrangian formulation for arbitrary complex fluids. 
Although the choice of the phase space T is in principle arbitrary, in the case of incom- 
pressible isothermal fluids, it is found [3| that the phase-space T can always be reduced to 
the direct-product space T = Q x V (restricted phase-space), where Q, V C IR 3 , Q is an 
open set denoted as configuration space of the fluid (fluid domain) and V is the velocity 
space. This type of approach (based on a restricted phase-space IKT formulation) will be 
adopted also in the sequel. The main motivation [of this work] is that some of the gen- 
eral understanding recently achieved in simple flows by means of the IKT approach could 
also give a significant contribution to a wider range of problems. In the sequel, we will 



concentrate on the issue of a consistent Lagrangian formulation for fluid dynamics based 
on a phase-space (IKT) description of incompressible fluids, whereby its pressure, velocity 
(and possibly also thermal) fluctuations are consistently taken into account in the phase- 
space dynamics of suitable phase-space inertial particles, i.e., particles whose dynamics is 
determined by the phase-space Lagrangian characteristics. The motion of these particles, 
as hey are pushed along erratic trajectories by fluctuations of the fluid-field gradients (in 
particular, characterizing the fluid pressure and temperature), is fundamental to transport 
and mixing processes in fluids. It is well known that the interaction between (deterministic 
and/or turbulent) fluctuations of the fluid fields and these particles still escapes a consis- 
tent theoretical description. Being a subject of major importance for many environmental, 
geophysical and industrial applications, the issue deserves a careful investigation. A key 
aspect of fluid dynamics is the correct definition of the (phase-space) Lagrangian dynamics 
which characterizes incompressible fluids. The customary approach to the Lagrangian for- 
mulation is based typically on a configuration-space description, i.e., on the introduction of 
the configuration-space Lagrangian characteristics r(t), spanning the fluid domain Q. Here 
r(t) denotes the solution of the initial-value problem: 



S = V(r,«), 



r(£ ) = r„, 

with r an arbitrary vector belonging to the closure £1 of £1 and V(r, t) being the velocity 
fluid field, to be assumed continuous in Q and suitably smooth in Q. In a previous paper 
[II] IKT has been proven to advance in time self-consistently the fluid fields, i.e., in such a 
way that they satisfy identically the requires set of fluid equations. For isothermal fluids, 
this conclusion is consistent with the results indicated elsewhere by Tessarotto et al.[A]. In 
particular, basic feature of the IKT approach is to permit to advance in time the fluid fields 
only by means of suitable evolution equations, without requiring the (numerical) solution 
of the Poisson equation for the fluid pressure. The purpose of this paper is to point out 
that, based the same approach [l|) for incompressible thermofluids, the dynamics of tracer 
particles, of finite mass (tut), which are injected in the fluid with arbitrary velocity, can be 
rigorously established. 



III. THE EXACT FRICTIONLESS CONDITIONAL DYNAMICS OF TRACER 
PARTICLES 

In this Section we intend to show that the equations of motion for a tracer particle in 
an incompressible fluid can be uniquely specified, once the velocity probability density are 
prescribed (conditional dynamics) . For definiteness in the sequel we shall consider an incom- 
pressible thermofTuid, described by the fluid fields {p = p > 0, V,p > 0, T > 0} , denoting 
respectively the (constant) mass density, fluid velocity, pressure and temperature describing 
the fluid. The fluid fields are assumed to be suitably smooth, classical solutions of the in- 

n 

compressible Navier-Stokes- Fourier equations (INSFE problem, |1|), all defined in a bounded 
connected set (fluid domain) 1]CR 3 (with f2 denoting its closure set where the mass density 
p is strictly positive): 

V • V = 0, (4) 

D - V = -— [Vp - f] + z/V 2 V, (5) 
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Dt 2c p \dx k dxij p c p 

In standard notation, -^V denotes the fluid acceleration, with -^ = J^+V- V the convective 
derivative, J is the quantity of heat generated by external sources per unit volume and unit 
time and finally f is the force density for which the Boussinesq approximation is invoked. 
Hence, f is taken on the form f = p g (I — k p T) + fi, where the first term represents the 
(temperature-dependent) gravitational force density and the second one (fi) the action of 
a possible non-gravitational externally-produced force density. As a consequence, in such 
a case the force density f reads f = p g (1 — k p T) + f 1( where the first term represents the 
(temperature-dependent) gravitational force density, while the second one (fi) the action 
of a possible non-gravitational externally-produced force. Hence f can be written also as 
f = — V0 + f.R, where = p gz and f^ = — p gk p T + fi denote respectively the gravitational 
potential (hydrostatic pressure) and the non-potential force density. Finally, u,x, c p and 
k p are all real positive constants which denote, respectively, the kinematic viscosity, the 
thermometric conductivity, the specific heat at constant pressure, and the density thermal- 
dilatation coefficient. A remarkable aspect of fluid dynamics is related to the construction 
of IKT's for hydrodynamic equations in which the fluid fields are identified with suitable 
moments of an appropriate kinetic probability distribution. This is achieved introducing a 



phase-space classical dynamical system (CDS) 



x D — ► x(£) — T tjio x , 



(7) 



which uniquely advances in time the fluid fields by means of an appropriate evolution oper- 
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Mo 



The CDS is assumed to be generated by the initial-value problem associate 



to the Lagrangian equations 



!r(i)=v(t), 
|v(t) = F(r(i),t;/), 



r(«o) 

V(to) 



(8) 



where x = (r, v) G Y and respectively r and v denote an appropriate state-vector and 
the corresponding configuration and velocity vectors. Moreover X(x, t) = {v, F(x,£;/)} 
is a suitably smooth vector fields, with F(x, t; f) representing an acceleration field. It is 
assumed that generally F(x, £;/) can depend functionally on f = /(x,£), to be identified 
with a suitable probability density function (pdf). In particular, let us require that x spans 
the phase-space T = QxV, where V = M. 3 denotes the velocity space. Therefore, introducing 
the corresponding velocity pdf /(x,t), defined so that 



d 3 vf(x,t) = 1, 



(9) 



v 



in T it fulfills necessarily the integral Liouville equation 



J(t)f(x(t),t) = f(x ,t ) 



(10) 



Here, x(t) denotes the Lagrangian path determined by the CDS ([Z|), which is assumed to 
be a diffeomeorphism at least of class C^ 2 '(T x I x I), with Id denoting an appropriate 
finite time interval, and J(t) 



dx(t) 
9x Q 



is corresponding Jacobian. Then, if the initial pdf 
/(x , t ) is suitably smooth, it follows that the time-evolved pdf /(x,t) satisfies necessarily 
the differential Liouville equation 

L/(x,t) = 0, (11) 

where L denotes the Liouville streaming operator. This equation (inverse kinetic equation), 
which may be interpreted as a Vlasov-type kinetic equation, can in principle be defined 
in such a way to satisfy appropriate constraint equations. In particular, thanks to the 
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arbitrariness of the dynamical systems (JTj), the velocity moments of /(x,£) can be identified so 
that suitable velocity moments of / coincide with the relevant fluid fields which characterize 
an incompressible thermofluid. In particular, introducing the kinetic pressure p\, defined so 
that 

Pi=Po(t)+p-<fi + poT/m, (12) 

we can identify m with the average molecular mass of the classical molecules forming the fluid 
and the ratio p /m = n with the fluid number density n of the fluid. In such a case, im- 
posing the constraints Y,pi— J d 3 vCr(x,£)/(x,£), respectively for G(x,t) = v,p D (v — V) 2 /3, 

v 
and subject to the requirement of strict positivity and regularity for /(x,t), the correct form 

of the acceleration field F = F(x, t; f) can be proven to be [l|: 

F(x,t;/) = F + F 1 , (13) 

where F , F x read respectively 

F (x,t; f) = - [V-S - V Pl + f R ] + u • VV+z/V 2 V, (14) 

Po 

F 1 (x,t;/) = i u (lA+lv-Q-l[V-a-Q} + ^V-n{^-|}, (15) 
2 [ Pl Pl p{ — ' J 2pi - [< 2 J 



and 

£>_ 1 



^V-— fp-z/V 2 V 



£>* Po 



+ ^, (16) 

m 
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dv, dv k Y i 



d 
A--^( Po +p)-V- 



2c p \da; fc dijl p c p 
Here the notation is given in accordance with Ref. [l]. Thus, po(t) (to be denoted as 
pseudo-pressure [4j) is a suitably-defined, strictly positive and smooth real function, to be 
determined to assure the validity of an H-theorem (and hence the strict positivity of /(x,t), 
Q and II denote the additional velocity moments Q = p J d 3 vu^-f and II = p J d 3 vuuf, 

The following consequences are implied by the previous Lagrangian equations (G|) [together 
with the definitions [Hfy - flty) ]: 

1. the Lagrangian equations (jHJ) can be interpreted either as deterministic or stochastic, 



depending whether the fluid fields themselves are treated as such (see also Ref. [261]); 



2. F(x, £;/) can be interpreted as the "conditional" Lagrangian acceleration, which de- 
pends functionally on the form of the pdf /(x,t). In particular, in the case of a 



turbulent fluid also the pdf must be considered as stocastic [26] 



3. it is possible to prove that the functional form of the vector field F(x, £;/) can be 
uniquely specified once /(x,t) is prescribed; 

4. the form of /(x,t) is determined uniquely by its initial condition /(x , t ) [see Eq. flTUl) ] 
and by the time-evolution operator T ijio ; 

5. in turn, T t)to is defined uniquely by the CDC (j7j) /(x,t); 

6. a particular solution of Eq. ffTTj) is delivered by the Maxwellian pdf 



/ M (x,t) =■ exp<^ — ^—}, ( 18 ) 



n v th,p I <4,p 



where u = v — V(r, £) and u^, p = ^/2pi(r, t)/p are respectively the relative and the 
thermal velocities. 

7. the case of an isothermal fluid is recovered by imposing that there results identically 
T(r,t) = const, and K = [see the r.h.s. of Eq.(JB])]. 

We conclude that, the Lagrangian equations |2P, when letting m = mr, define the 
equations of motion for a tracer particle with initial conditions {r(t ), v(t )} = (r , v ) and 
subject to the conditional acceleration F(x, t; f). 

IV. CONCLUSION AND OUTLOOK 

In this Note we have proven that the exact dynamics of tracer particles injected in in- 
compressible thermofluids can be rigorously established. The result appears relevant not 
only from the theoretical viewpoint but also for its applications and in particular for the 
description and monitoring of pollutant dynamics in the atmosphere. 

We stress that only weak restrictions have been assumed on the form of the pdf, affecting 
the form of the conditional acceleration F(x, t; /). In addition, the theory applies in principle 
to arbitrary fluid fields which are strong solutions of the INSFE problem. Therefore, this 
permits in principle also the proper treatment of turbulent, i.e., stochastic fields. As a 

10 



consequence, the formulation here presented is susceptible of applications in turbulence 
theory. 
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